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ABSTRACT 

We study statistical properties of 3D classical spin glass layer of certain width and 
infinite length. The 3D spin glass is represented as an ensemble of disordered ID 
spatial spin-chains (SSC) where interactions are random between spin-chains (non- 
i-^ I ideal ensemble of ID SSCs). It is proved that at the limit of Birkhoff's ergodic 

1) ' hypothesis performance 3-D spin glasses can be generated by Hamiltonian of disor- 

dered ID SSC with random environment. Disordered ID SSC is defined on a regular 
lattice where one randomly oriented spin is put on each node of lattice. Also it is 
supposed that each spin randomly interacts with six nearest-neighboring spins (two 
spins on lattice and four in the environment). The recurrent transcendental equa- 
te ' tions are obtained on the nodes of spin-chain lattice. These equations combined with 
the Silvester conditions allow step by step construct spin-chain in the ground state 
'^ ■ of energy where all spins are in minimal energy of classical Hamiltonian. On the 
^ , basis of these equations an original high-performance parallel algorithm is developed 
Q I for 3-D spin glasses simulation. Distributions of different parameters of unperturbed 
spin glass are calculated. In particular, it is analytically proved and by numerical cal- 
culations shown that the distribution of spin-spin interaction constant in Heisenberg 
nearest-neighboring Hamiltonian model as opposed to widely used Gauss-Edwards- 
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have variance. A new formula is proposed for construction of partition function in 
kind of one- dimensional integral on energy distribution of 1-D SSCs. 
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I. INTRODUCTION 

The wide class of phenomena and structures in physics, chemistry, material science, biol- 
ogy, nanoscience, evolution, organization dynamics, hard-optimization, environmental and 
social structures, human logic systems, financial mathematics, etc. are mathematically well 
described in the framework of spin glass models [l|-l9J. 

The considered mean-field models of spin glasses as a rule are divided into two types. The 
first consists of the true random-bond models where the coupling between interacting spins 
are taken to be independent random variables 1CH12I] . The solution of these models is ob- 



tained by n-replica tric k |10l . |12| and requires invention of sophisticated schemes of replica- 
symmetry breaking [12, ll3| . In the models of second type the bond-randomness is expressed 
in terms of some underlining hidden site-randomness and is thus of a superficial nature. 
It has been pointed out in the works |l4|-|l6|. however, this feature retains an important 
physical aspect of true spin glasses, viz. they are random with respect to the positions of 
magnetic impurities. 

Note that all mentioned investigations as a rule conduct at equilibrium's conditions of 
medium. This fact plays a key role both in analytical and numerical simulation by Monte 
Carlo method. 

Recently, as authors have shown [l7| some type of dielectrics can be studied by model of 
quantum 3-D spin glass. In particular, it was proved that the initial 3D quantum problem 
on scales of space-time periods of an external fields can be reduced to two conditionally 
separable ID problems where one of them describes an ensemble of disordered ID spatial 
spin-chains between which are random interactions (further will be called nonideal en- 
semble). 

In this paper we discuss in detail statistical properties of classical 3D spin glass with sug- 
gestion that interactions between spins have short-range character. We prove that nonideal 
ensemble of ID SSCs exactly describes the statistical properties of classical 3D spin glasses 
in the limit of Birkhoff's ergodic hypothesis performance. In the work a new high perfor- 
mance algorithm for simulation of this traditionally difficult calculated problem is developed. 

In section II the classical spin glass problem on 3D lattice is formulated. Equations 
for stationary points and corresponding Silvester conditions are obtained for definition of 
energy minimum on lattice nodes (local minimum of energy). The formula for computation 
of different parameters distributions of spin glass is defined. 

In section III a theorem on reduction of 3D spin glass problem to the problem of nonideal 
ensemble of ID SSCs is proved. 

In section IV numerical experiments are adduced for unperturbed ID SSCs ensemble with 
spin-chain's length 10^ (Iq. In particular, distributions of energy, polarization and spin-spin 
interaction constants of nonideal ensemble are investigated in detail. 

In section V partition function is investigated in detail in the configuration integral's 
representation. A new representation is suggested for partition function in kind of one 
dimensional integral on energy distribution of nonideal ensemble. 

In section VI the obtained theoretical and computational results are analyzed. It is very 
important to note that it has been proved that in the framework of the developed method 
it is always possible to exactly compute the ground state energy of 3D spin glasses. 








■♦ X 



a P 

FIG. 1: The structure of amorphous sihcon dioxide Si02 is described by 3-D random network with 
covalent bonds. Every sihcon vertex (gold sphere) has 4 edges and every oxygen vertex (red sphere) 
has 2 edges. 

II. FORMULATION OF PROBLEM 

The objects of our investigation are solid-state dielectrics, type of Si02 glass (amorphous 
silicon dioxide). According to the numerical ab initio simulations |7|, the structure of this 
type compound can be well described by 3D random network (Fig. la). The red and brown 
lattice points on this figure correspond to different atoms, while the links between them 
correspond to covalent bounds. As a result of charges redistribution in outer electronic 
shells, atoms of Si acquire the positive charge and atoms of O correspondingly the negative 
charge. Thus, we can consider compounds of this type as a disordered 3D system of similar 
rigid dipoles (hereinafter termed as a system of 3D disordered spins. Fig. lb). Let us remind 
that under the similar rigid dipoles are meant the dipoles for which the absolute values are 
equal (IpJ = |pj| = p^, where pi and pj are two arbitrary dipoles), and they don't vary 
under the influence of an external field. 

The Hamiltonian of 3D classical spin glass system reads: 



^({r}) 



7 ^ JijSiSj, {r} = ri,r2; 
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where indices i and j run over all nodes of 3D lattice, r^ correspondingly denotes the co- 
ordinates of 2-th spin (see Fig. lb). For further investigation we will consider a spin glass 
layer of certain width L^ and infinite length (see Fig. 2). We will consider 3D compound 
in the framework of nearest-neighboring Hamiltonian model. Let us note that even for this 
relatively simplest model numerical simulations of spin glasses are extremely hard to solve 
NP problems. 

At first we will consider an auxiliary Heisenberg Hamiltonian of the form: 
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FIG. 2: ID SSC with the random environment. Recall that each spin-chain is surrounded by four 
spin-chains which are randomly interacted with it. By symbols ® are designated spins from the 
random environment (four spin-chains of surrounding) . 
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where the first term Hq ({r}; A^, 
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H^'\{r}; N,) = - J^ J^^+lS,S,^u 



i=0 



describes the disordered ID spatial spins chain (SSC) while the second term Hq ({r}; N. 
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correspondingly describes the random surroundings of ID SSC (see Fig. 2). In ([T]) Ja+i and 
Jii^ are correspondingly random interaction constants between arbitrary i and z + 1 spins 
and between i and i^ spins, Si, Si+i and 5j^ are spins (vectors) of unit length, which are 
randomly orientated in 0(3) space. From the general reasons it follows that with the help 
of ([T]) Hamiltonian and by way of successive constructing we can restore the Hamiltonian 
of 3-D problem. Recall that the meaning of the construction is as follows. On the first 
step the central spin-chain on the x-axis with its surroundings from four random spin- 
chains is considered (see Fig. 2). On the second step as central spin-chains are considered 
corresponding spin-chains from the random surroundings each of which are surrounded by 
new four neighboring spin-chains. Thus, repeating this cycle periodically we can construct 
the Hamiltonian of 3D problem. This idea will be rigorously proved below. 

For further investigation of spin glass problem it is useful to write Hamiltonian ([T]) in 



spherical coordinates system: 



i=0 



Ho{{r}; N^) = - ^S Ja+Acos xpi cos i/ji+i cos{(pi - (fii+i) + sin ^j sin ^i+i 

4 

y^ Jii^ cos ^i cos iJi^ cos{ipi -ifi^) + sin ^pi sin i^i 



a=\ 



(2) 



Now the main problem is to find the angular configurations and spin-spin interaction con- 
stants which can make the Hamiltonian minimal on each node of lattice. 
Let us consider the equations of stationary point: 






0, 






0, 



(3) 



where 0j = {ipi, ipi) defines the orientation of i-th spin {ipi, ipi are correspondingly the polar 
and the azimuthal angles). In addition, = (6i, 62....6Ar^) describes the angular configu- 
ration of spin-chain consisting of N^ spins. 

Substituting ([2]) into ([3]) we can find the following recurrent equations: 
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sin^j cos^j_i cos{ipi — v^i-i) + cos^j sin-j/ij. 



Jii+i -sin^jCos^i+icos(v9j - ^i+i) + cos ^j sin ?/^i+i 



^Jii„ - sin i^i cos i^i^ cos{(pi - fi^) + cos i^i sin ^j^ 



+ 
= 0, 



(T=l 



Ji_ii cosV'j-i sin(v7i - ipi-i) + Jii+i cos V'i+i x 
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sin(v9i - (/?i+i)J + ^ Jiv cosipi^ sm{ipi - ipi^)j costpi = 0. 



(4) 
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In order to satisfy the conditions of local minimum (Silvester conditions) for Hq, it is 
necessary that the following inequalities are carried out: 

A^^^M) > 0, A^,^,(e°)A^,^,(e°) - 4^^^ (6°) > 0, (5) 

where A^^a, = d'^Ho/daf and A^^/b^ = A^^p^ = d'^Ho/daidf3i, in addition: 

^^«V'i(Qi') = ^i-iijcos^/^^cos^/'i-i cos((y9° - (fi-i) + sin ?/^° sin ^i_i I + Jii+i|cos^°cos?/'i+iX 
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cos(v3- - <fi+i) + sin ?/^° sin ^pi+i j + ^ -^ii, jcos ^° cos ^pt^ cos(v?° - fij + sin ^|)^ sin ^pi^ | , 

(T=l 

4 
\^V^^{'^'i) = \ji-liCOSl/ji-iCOs{ip°i - ifii-i) + Jii+iCOSIpi+iCOsiip'^ - ipi+l) + y^ JiyX 
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cosijj,^ cos((/?° -iPij] cos^°, A^,^,(e°) = 0. 



Recall that 6° = {ipi: V^?) designates the angular configuration of the spin in case when the 
condition of local minimum for Hq is satisfied. 

Thus, it is obvious that the classical 3D spin glass system (see Fig. lb) can be considered 
as an nonideal ensemble of ID SSCs (see Fig. 2) and there are random interactions between 
spin- chains. 

Now we can construct distribution functions of different parameters of ID SSCs nonideal 
ensemble. To this effect it is useful to divide the nondimensional energy axis e = e/6e into 
regions > Eq > ... > En, where n >> 1 and e is the real energy axis. The number of 
stable ID SSC configurations with length L^ in the range of energy [e — 6e,e + Se] will be 
denoted by Ml^{e) while the number of all stable ID SSC configurations - correspondingly 
by symbol M[^ = J21=i^L^i^j)- Accordingly, the energy distribution function can be 
defined by the expression: 

FL.Ae;do{T)) = MUe)/M[:\ (6) 

where distribution function is normalized to unit: 

n „0 

lim ^Fz.,(£,;c?o(T))fe, = / ^^^(e; c?o(T))& = 1. 

By similar way we can construct also distribution functions for polarizations, spin-spin 
interaction constant, etc. 



III. REDUCTION OF 3D SPIN GLASS PROBLEM TO ID SSCS ENSEMBLE 

PROBLEM 

Modeling of "iD spin glasses is a typical NP hard problem. This type of problems are 
hard-to-solve even on modern supercomputers if the number of spins in the system are 
more or less significant. In connection with told, the significance of new mathematical 
approaches development is obvious and on the basis of which an effective parallel algorithms 
for numerical simulation of spin glasses can be elaborated. 

Theorem: The classical 3D spin glass problem at the limit of isotropy and homogeneity 
(ergodicity) of superspins distribution (sum of spins in chain) in 3D configuration space is 
equivalent to the problem of disordered ID SSCs ensemble. 

It is obvious that the theorem will be proved if we can prove that in case when the distri- 
bution of superspins in 3-D configuration space is homogeneous and isotropic, the following 
two propositions take place: 

a) In any random environment which consists of four arbitrary spin-chains it is always pos- 
sible to find at least one physically admissible solution for spin-chain (the direct problem), 
and 

b) It is possible to surround an arbitrary spin-chain from the given environment with such 
environment which can make it physically admissible spin-chain solution (the reverse prob- 
lem). 

The direct Problem. 
Using the following notation: 

^i+i = cos tpi+i, rii+i = sin(v?i - (Pi+i), (7) 



equations system ([6]) can be transformed to the following form: 



Ci + Jh+i [y 1 - if+i - tan^i ii+l^Jl - vf+i] = 0' 

where parameters Ci and C2 are defined by expressions: 

4 
Ci = Ji^ii[smipi_i - tan?/^iCosV'i-i cos{ipi - (fi^i)] + ^ Jj 

sinT/ij^ - tan^/'j cos^j^ cos(v?i - (fij 



(8) 



o-=l 



C2 = Ji-iiCOSipi_i sm{ipi - ipi_i) + ^ Jiv cos^i^ sm{ipi - ipi^ 



a=l 



From the system (jS]) we can find the equation for the unknown variable r^j+i: 



Cir/,+1 + C2JI- vf+i tanipi + J Jf^^^rit+i - C^ = 0. 



(9) 



We have transformed the equation Qj to the equation of fourth order which is exactly solved 

further: 
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(10) 



where 



A = Cl[ 4+1 cos^ ^i + C^ + 2Cl sin^ ^, [l ± C^' y^Jfi+i - Cl - C| cot ^,] } . 






C3 = -Ci + C2' sin' z^„ E = 4+1 cos^ ^i + 2C3 J;,+i cos' ^^ + (Cf + C2' sin' 7/;, 

Note that from the condition of nonnegativity of the value under the root we can find the 
following nonequality: 



4m > ^1+^2- 
In consideration of (ITj), we can write following conditions: 

< Cli < 1, < r,^^ < 1. 



(11) 



As it follows from equations ( jTOj) . if the solutions in previous two nodes {i — 1) and i are 
known, then the solutions {ipi+i,(pi^i) in the node {i + 1) can be defined only by constant 
Jii+i- In this connection a natural question arises - are there solutions for spin-chain in 
arbitrarily given environment? 

Let us consider Silvester conditions (E]) which can be written in the form of the following 
inequalities: 



Jii+i cos^°cos^i+i cos(v3° - (Pi+i) > -fli - sin^°sin?/;i+i, 
Jii+i cosV^j+i cos(v9- - (pi+i) cosipi > -as, 

where constants Oi and 02 are defined by expressions: 

4 
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cos?/;°cos'?/'j_icos(v9° - v9j_i) + sin?/;°sin'?/'j_i +22'^"'^ cos^°x 
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FIG. 3: The projection of spin-chains ensemble on the {Y, Z) plane. Spin-chains are designated by 
symbols Q ^'^d (^ which correspondingly form the old and new environments. 



cosil^i^ cos(v2° - </?v) + sin ?/;° sin ^j, 



4 

02 = |ji_iiCOS?/^i_iCOs(v?- -</?i_l) + ^ Jii^ COS Ipi^ COs{ip^i -ipi^)j COS ij'^. 



a=l 



So, the problem leads to the answer of the following question - are inequality (TTT]) and (TT2|l 

compatible or not. Taking into account solutions (fTO!) it is easy to prove that conditions (TT2|) 

are automatically compatible at large absolute values of Jjj+i. On the other hand, there is 

no any contradiction with condition fill I) . Thus the direct problem or the proposition a) is 

proved. 

Now our aim is to prove the reverse problem or the proposition b) which consists of the 

following. We choose a spin-chain from the environment (see Fig. 2), for example {^4} = 



(4) 



.(4)^ 



(O4, I4, ..., Nxi). In this spin-chain all angular configurations of spins (69 , ...6]v^) are known 
but the constants that define spin-spin interactions in spin-chain and interactions between 
spin-chain and its environment still are not defined. We will prove that it is always possible 
to surround each spin-chain by such environment that the selected spin-chain will be the 
correct solution from the main physical laws point of view (see conditions (j4])-([5])). In 
the considered case {i^} = {i^} spin-chain is surrounded by four neighbors, one of which 
{io} = {^2} is fully determined while three spin-chains {«i},{«3} and {i^} should be still 
specified (see Fig. 3). Recall that the mark " " designates a new environment with three 
spin-chains. However, for simplicity we will omit or more clearly make change them in the 
subsequent calculations {{io}, {h}, {h}^ ih}^ ih}) ~^ ({^o}, {^i}{^2}, {is}, {u})- The proof 
of the proposition should be conducted as follows. We will suppose that the constants of 
spin-spin interactions in considered chain and corresponding parameters of two spin-chains 
of environment are known. We will show that by special choosing of parameters of the third 
spin-chain {^3} it is possible to ensure the condition of local minimum energy is satisfied in 
the considered spin-chain. 



Let us define the following denotations for constants: 

ci = Ji-ii[- sin ipi cos ipi-icos{(fi - v?j-i) + cos ipi sin ipi^i] + 

Jii^ [-sin^iCos^i^cos(v9i - fij + cos^j sinz/;^^], 

C2 = - sin ipi cos ^i+i cos(v3i - ^i+i) + cos ipi sin -j/^j+i, 

C3 = Ji__iiCos^i_isin(v9i -<fi-i) + Jii^ cos tpi^ sin{Lpi - Lpi^), 

C4 = cosipi+i sin{(fi - (fi+i), a = A. (13) 

Using flT3l) we can transform equations (jl]) to the following form: 

3 

ci + C2Jii+i + '^ Jii^[- sini}iCosi)i^ cos{ipi - (fi,) + cos^i sin^j^] = 0, 

(T=l 

3 

C3 + C4Jii+i + ^Jii„ COS ipi^ sin(v?i - v^v) = 0, 

(T=l 

which are equivalent to the following relations: 

1 ^ 
Jii+i = — y^ Jii^[-sin?/;iCos^i^cos(v9i- v^i,) + cos ipi sin ip^^], 



3 

C4 C4 



1 
Jii+i = — ^ Jii, costpi^ sin{ipi - ipi^). (14) 



cr=l 

After excluding Jjj+i from (IT^ we find the following equation: 
3 



V's -^[-sinV'jCosV'i, cos(v2i-(^i^) + cos V'i sin V'i J - 

1 I C2 

a=l ^ 

"" cosilJi^sin{(pi-ipi^)} -C5 = 0, C5 = — - — . (15) 



C4 J C2 C4 

Having made the following designation: 

2 



_D = y^J ^!f.[_sin^-cos^i^cos(v2i - (fi^) + cos tpi sin ipi^] 

1 I C2 

cosipi^sin{ipi-ipi^) > -C5, 



cr=l 

C4 

we can transform equation (TTSl) to the following form: 

-D H — ^[— sin '?/'j cos -i/^ig cos{ipi — ipi^) + cos-^/'jSin-i/^ig] 
C2 

^ cos V'ia sin{(pi -fi.,) = 0. (16) 

C4 

Now substituting: 



X = cos ^jJi^, (17) 



10 

in ( !T6|) we find the equation: 

D H — ^[— xsin^j cos{ipi — (pi^) + Vl — x"^ cos^j] 

-a;^^sin(v9i -(/Pig) =0, (18) 

C4 

From (ITSl) the following square equation can be found : 

Kox^ + 2Kix + i^2 = 0, (19) 

where the following designations are made: 



Kq = cos^ ipi + I sin^j cos{ipi - (pi.^) -\ sm{ipi - (p. 



2 

siniipi - ipi^) 1 
C4 



Ki = — -^ ( sin ?/^i cos((y9i - v^ig) + —sm{ipi - ^i^) ), K2 = { —^ ) 
Discriminant of the square equation ( 1T9|) has the form: 



2 



cos V'i 



Da; = ( sin V^i cos(v2i - V^ig) + — sin(v9i - v^ig) j cos^ ^i 

+ |cos2 V-i - (^j I cos^^Ai > 0, (20) 

which on some set of Ja^ can be positive, i.e. z-th spin in spin-chain {^4} will satisfy the 
local minimum conditions. 
Let us define: 

y = cos{^i-^i,,), (21) 

Substituting ( 1211) in ( IT6l) we will find that: 

D H — ^— [—y sin ipi cos ipi^ + cos ipi sin -j/^jg] ''— cos -j/^jg v^l — y'^ = 0, 

C2 ' C4 

After squaring we will have the following equation: 

Moy^ + 2Miy + M2 = 0, (22) 

where the following designations are made: 

C2V . . 2 l\ 2 . n. . . , ( , . , , ^C2 



Mo =11 — ) + sin '?/'j ) cos -j/^jg , Mi = — sin t/^j cos ■i/'jg } cos ■i/'j sin -i/^jg + 



. , rti --"- rtg, _L T t. rtg \ "-- rt r tg i y 



'"3 



-^C2 \ / C2 



M2 = 1 COS ipi sin ijji + -; — ) — 1 — ) cos ipi^ . 



2 

I (:( )s 
^i is, J \ '^4 , 



'lis 

The discriminant of the square equation ( l22l) has the form: 

Dy = { — j cos^ ifji + sin^ ipi cos^ ipi^ — I — h cost/'j sin-i/'jg J > 0. (23) 



11 

Obviously there are some set of constants Ja^ on which Dy > 0. However, it is more impor- 
tant to find the region of the interaction constant Jjjg values for which both determinants 
D^ and Dy are positive. 

In particular as the analysis of the following condition shows: 



Dco 



cos Ipi 



— ^iii — 



DC2 



COS Ipi 



(24) 



discriminant D^ is always nonnegative. From the other side: 



sin^,3^-- ^, (25) 

Jii^COSljJi 

which will assure that Dy discriminant is always nonnegative. A simple analysis of conditions 
(IMl) and ( l25l) shows that they are compatible. In other words the set of constants Ja^ which 
satisfies the energy local minimum condition is not empty and therefore the proposition b) 
is proved. 

So, we have proved the validity of a) and b) propositions. It is obvious that at the 
simulation of ID SSC problem we can by this way fill up 3D space by ID SSC which is 
equivalent to obtaining 3-D spin glass. In case when the number of ID SSCs is so much 
that the directions of spins in 3-D space are distributed isotropically and homogeneous, the 
statistical properties of both problems (3-D spin glass and 1-D SSCs nonideal ensemble) will 
be obviously identic. 

The theorem is proved. 



IV. RESULTS OF PARALLEL SIMULATIONS 

One important consequence of the theorem is that for the numerical simulation of the 
problem we can use the algorithm for solving the direct problem. Obviously, a large number 
of independent computations of ID SSC which can be carried out in parallel and in statistical 
sense make it equivalent to the problem of 3-D spin-glass. This approach considerably reduces 
the amount of needed computations and helps us effortlessly simulate statistical parameters 
of 3-D spin glasses of large size. 

The strategy of simulation consists of the following steps (see Fig. 4). At first, the angular 
configurations of four spin-chains are randomly generated which form random environment of 
the spin-chain which we plan to construct later. On a following step a set of random constants 
Jii^ are generated, which characterizes the interactions between the random environment 
and the spin-chain. The interaction constants are generated by Log-normal distribution. 
The angular configurations of the random environment are generated the same way as it is 



described in [18[. Now when the environment and its infiuence on disordered 1-D SSC are 
defined, we can go over to the computation of spin-chain which must satisfy the condition of 
local energy minimum. Note that the scheme of further computation of nonideal ensemble 
of ID SSCs (see Fig. 2) is identical to the scheme of the computation of an ideal ensemble 



of disordered 1-D SSCs (see [l8|). Note that all calculations of 1-D SSCs nonideal ensemble 
are done for spin-chains with 10^ do length which require huge computational resources. 

As the simulations show, for the ensemble which consists of 10^ spin-chains, the dimen- 
sional effects practically disappear (see Figs. 5a, 5 & and 6) and the energy distribution F{e) 
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FIG. 4: The algorithm of parallel simulation of statistical parameters of disordered ID SSCs 
nonideal ensemble. The symbol fJ^ describes the input of environment, M is a number of simulation 
or overall number of spin-chains in the nonideal ensemble, N^ is a number of spins in chain. 
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FIG. 5: On figure a) the energy distribution of \D nonideal ensemble of SSCs with 10'^ length is 
shown. The red line shows a numerical data visualization while the green one illustrates its fitting 
by Gaussian function. On figure V) the visualization of numerical data of spin-spin interaction 
constants (pink line) is shown while the blue one denotes Gaussian distribution. The analysis of 
the numerical data proves, that the green curve is not analytic function and by the character is 
the Lev's skew a-stabile distribution function. 
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FIG. 6: On figure the polarization distributions on different coordinates are shown after 10^ sim- 
ulations. 



has one global maximum and is precisely approximated by Gaussian distribution (see Fig. 
5a). 

Mean values of polarizations on coordinates are not very small, especially when it comes to 
coordinate x (thickness of spin glass layer): p^ = —0.13508, Py = 0.036586, Pz = —0.059995 
and correspondingly the average energy of 3D SSC is equal to e = —990.88, where p = 
J_^ F{p)pdp, p = {px, Py, Pz), £ = /„oo F{e)ede and F is the distribution function. As our 
numerical investigations have shown on the example of systems where thickness of spin glass 
layer is not so large oc 25do -^ lOOdo, for a full self-averaging of superspin it is necessary to 
make oc N^ simulations. In other words, the system can be fully ergodic in considered case 
if we continue the numerical simulations of the spin-chains up to oc 10^ times. 

It is analytically proved and also the parallel simulation results show that the spin- 
spin interaction constant cannot be described by Gauss-Edwards- Anderson distribution (see 
Fig. 5b). It essentially differs from the normal Gaussian distribution model and can be 
approximated precisely by Levy skew alpha-stable distribution function. Let us recall that 
Levy skew alpha-stable distribution is a continuous probability and a limit of certain random 
process X{a, /3, 7, 6; k), where the parameters correspondingly describe an index of stability 
or characteristic exponent a G (0;2], a skewness parameter /3 G [—1; 1], a scale parameter 
7 > 0, a location parameter 5 G R and an integer k which shows the certain parametrization 



see 



19, Ml). Let us note, that the mean of distribution and its variance are infinite. 



However, taking into account that spin-spin interaction constant has limited value in real 
physical systems, it is possible to calculate distribution mean and its variance. In particular 
if J G [-5, +5] then J = 0.89717 and 1^ = 5.3382. 

In the work are also presented polarization distributions on different coordinates (see Fig. 
6). As for the polarization distributions, they are obviously very symmetric by coordinates 
in the considered case (see Fig. 6). 

One of the advantages of the developed algorithm is that we are able to take into account 
the branching solutions at the successive constructing of the spin-chain (see Fig. 7). As 
calculations show, the number of branching solutions z/ for spin-chains of length 10^ do is 
not more than 25. At the simulation process only those spin-chains are considered for 
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FIG. 7: On figure the number of brandling of solutions v is shown along with the spin-chain 
length depending on different initial conditions which are indicated by way of various colors. 

which Silvester conditions are satisfied on each node. If on some node the conditions are not 
satisfied we try to regenerate Jjj+i in order to obtain a new solution. However, if the solution 
is not found after large quantity simulations it means that the weight of these solutions all 
are extremely small and further simulations of these spin-chains are unpractical. Thus, 
in case when the ensemble consists of a large number of spin-chains, the self-averaging of 
superspin (sum vector of spin-chain) in 3-D space occurs with high accuracy. It is important 
to note, that the summation procedure on the number of spins in chain or on the number of 
spin-chains in ensemble, is similar to the procedure of averaging by the natural parameter 
or "timing" in the dynamical system. The latter means that at defined space scales of spin 
glasses it is possible to introduce the concept of ergodicity for both separate spin-chains and 
ensemble as a whole. 



V. PARTITION FUNCTION 



The main object of investigation of statistical mechanics, information science, probability 
theory and etc., is the partition function which is defined for classical many particle case in 



configuration space as follows [21 



Z{fi) 



exp[— /3iJ({r})]dri(ir2... 



/S 



knT 



(26) 



where /cs is the Boltzmann constant and T is the thermodynamic temperature. Obviously, 
when the number of spins or spin-chains in the system are large we can consider integral 
( l26l) as a functional integral. In any case the number of integration in expression ( l26ll as 
a rule is very large for many tasks and the main problem lies in the correct calculation 
of this integral. However, in the representation of f l26|) configurations of spin-chains that 
are not physically realizable obviously make a contribution. Moreover, the weight of these 
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configurations is not known in general scenario and it is unclear how to define it. With this 
in mind and also taking into account the ergodicity of the spin glass in the above mentioned 
sense, we can define the partition function as: 



Z,{(3-N, 



exp \^j3e\ F{e;Nx)de, 



(27) 



where F{e; N;j.) is the energy distribution function in nonideal ensemble of ID SSCs with 
certain length A^^, (see also definition ([6])). 

Now we can define the Helmholtz free energy for ensemble of ID SSCs by two different 
ways. Using standard definition for Helmholtz free energy we can write: 



Q(/3;iVx 



iV./3 



\n[Z{/3;N^)], Q.(/3;iV, 



iVx/3 



ln[Z,(/3;iV,)]. 



(28) 



Note that the dependence on N^ of the expressions in (l28l) arises due to the finite layer 
width. In particular, using the expression of partition function (|26l) we can find the average 
value of free energy coming on one spin in the chain (see also |22l|): 



Q(/3;iV.) 



N,/3 



Na:-l 



j=0 



sinh Xj 



Xi 



Xj, 



Jii+lP 



(29) 



where (...) designates averaging by ID SSCs ensemble. Now the main problem is the 
investigation of behavior of free energy subject to the parameter (3. Correspondingly we can 
define two forms of free energy derivatives: 



g(/3;iVx) 



dp 



N.P' 



N^-l 



i=0 



E ! + '■> 



sinh Xj 



Xi 



Xi coth Xi 



q*{/3]N,) 



dQMN.) 



0(3 N^pyp 

where Z,.^(/3; N,) = dZ,{l3- N,)/dl3. 



' ^^nZ.(/5;iV.)-^-(^^^^ 



Z,(/3;iV,) 



(30) 



VI. CONCLUSION 



A new parallel algorithm is developed for the simulation of the classical "iD spin glasses. 
It is shown that 3D spin glasses can be investigated by the help of an auxiliary Heisen- 
berg Hamiltonian ([1]). The system of recurrent transcendental equations ([3]) and Silvester 
conditions (j4]) are obtained by using this Hamiltonian. Let us note that exactly similar 
equations of stationary points ([3]) also can be obtained if the full 3D Hamiltonian (see the 
first unnumbered formula) is used in the framework of short-range interaction model. That 
allows us step by step construct spin-chain of the specified length with taking into account 
the random surroundings. It is proved that at the limit of Birkhoff's ergodic hypothesis 
performance, 3D spin glass can be generated by Hamiltonian of disordered ID SSC with 
random environment. We have proved that it is always possible to construct spin-chain in 
any given random environment which will be in ground state energy (direct problem). We 
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FIG. 8: On the figure two derivatives of Helmholtz's free energy are shown depending on reverse 
temperature beta which have been calculated by formulas (j30p after 10® simulations. The figure 
shows that at low temperatures the curves approach each other which is rather natural, resulting 
from the decrease of the system's entropy. 



have also proved the inverse problem, namely, every spin-chain of the random environment 
can be surrounded by an environment so that it will be the solution in the ground state. In 
the work all the necessary numerical data were obtained by way of large number of paral- 
lel simulations of the auxiliary problem in order to construct all the statistical parameters 
of 3-D spin glass at the limit of ergodicity of ID SSCs nonideal ensemble. As numerical 
simulations show, the distributions of all statistical parameters become stable after oc A^^^ 
independent calculations which are realized in parallel. The idea of ID spin-chains parallel 
simulations, based on this simple and clear logic, greatly simplifies the calculations of 3D 
spin glasses which are still considered as a subset of difficult simulation problems. Let us 
note that computation of spin-spin interactions distribution function from the first princi- 
pals of the classical mechanics is very important result of this work. As analysis show, the 
distribution is not an analytic function. It is from the class of Levy functions which does 
not have variance J^ and mean value J. 

Despite the absence of calculations by other methods, it is obvious that the developed 
scheme of calculations should differ from other algorithms, including the algorithms which 
are based on Monte Carlo simulation method J23|, by the accuracy and efficiency. We 
were once again convinced in the accuracy and efficiency of the algorithm after analyzing 
the results of different numerical experiments by modeling the statistical parameters of 3D 
spin-glass system which are presented in figures 5 a, b and 6. 

In the work a new way of partition function construction (configuration integral) is pro- 
posed in the form of one- dimensional integral of the energy distribution, which unlike the 
usual definitions does not include physically unrealizable spin-chains configurations (see dif- 
ference of free energy derivatives on Fig. 8). It is obvious that the new definition of partition 
function is more correct and in addition it is very simple for computation. 

Finally, the developed method can be generalized for the cases of external fields which 
will allow us investigate a large number of dynamical problems including critical properties 
of 3D classical spin glasses. 
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